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We present the first calculation of coherent backscattering 
with inelastic scattering by saturated atoms. We consider the 
scattering of a quasi-monochromatic laser pulse by two distant 
atoms in free space. By restricting ourselves to scattering of 
two photons, we employ a perturbative approach, valid up to 
second order in the incident laser intensity. The backscatter- 
ing enhancement factor is found to be smaller than two (after 
excluding single scattering) , indicating a loss of coherence be- 
tween the doubly scattered light emitted by both atoms. Since 
the undetected photon carries information about the path of 
the detected photon, the coherence loss can be explained by 
a which-path argument, in analogy with a double-slit experi- 
ment. 

PACS numbers: 42.25.Fx, 32.80-t, 42.50-p 



I. INTRODUCTION 

Weak localization of light in random media was demon- 
strated for the first time in the eighties [1-3]. Here, con- 
structive interference between two waves which interact 
with the same particles, but in reversed order, enhances, 
in average, scattering in the direction opposite to the 
incident light. For systems obeying the reciprocity sym- 
metry [4], the backscattering enhancement factor, i.e. the 
light intensity detected in exact backscattering direction 
divided by the background intensity, is exactly two, pro- 
vided single scattering can be removed. For spherically 
symmetric scatterers, the latter is achieved in the helicity 
preserving polarization channel. 

Similar interference effects between multiply scattered 
waves also affect the properties of transport through dis- 
ordered media. If the mean free path can be sufficiently 
reduced, the transport is even expected to come to a com- 
plete standstill [5] . In experiments on strong localization 
of light [6], however, the role of absorption is discussed 
controversially [7,8]. 

One may wonder whether a medium consisting of in- 
dividual atoms would constitute a good candidate for 
strong localization. In contrast to the classical scenario 
(Maxwell's equations in a medium with random dielectric 
constant), the quantum- mechanical atom- photon inter- 
action exhibits some characteristic features, which may 
affect the coherence between multiply scattered waves. 
Firstly, the resonance may be extremely sharp, corre- 
sponding to a very narrow linewidth T of the excited 
state. On the one hand, this leads to a large atom-photon 
scattering cross section and slow diffusion of light [9] - 



properties in favor of localization. On the other hand, it 
implies that the atoms have to be cooled to very low tem- 
peratures. Only if the Doppler shift induced by a moving 
atom is much smaller than T, the interference between 
two counterpropagating waves is preserved [10]. Typi- 
cally, this regime is reached at about a few mK, which 
is, however, still high enough to neglect the thermal dc- 
Broglie wavelength of the atoms, i.e. to treat their exter- 
nal motion classically. Furthermore, atoms usually have 
an internal quantum structure, which may have a strong 
impact on coherent backscattering [11-14]. If necessary, 
this can be circumvented by using atoms with a nondc- 
generate ground state (J = 0) [15]. 

Another property of the atom-light interaction, whose 
impact on coherent backscattering has so far remained 
almost unexplored, is the strongly nonlinear response of 
an atom to incoming radiation. A single photon being 
sufficient to bring the atom to the excited state, where 
it rests for a quite long time T -1 without being able to 
scatter other photons, a saturation of the atomic medium 
can be induced already with rather moderate laser inten- 
sities. Not only the atom-photon cross section, but also 
the spectrum of the light is affected by saturation. With 
increasing saturation, it becomes more and more proba- 
ble that an atom scatters inelastically, i.e. that it emits 
photons at a frequency different from the one of the inci- 
dent laser. As we will show in this paper, this implies a 
loss of coherence between two reversed scattering paths. 
Similarly, a recent experiment showed coherent backscat- 
tering by a cloud of cold Strontium atoms to be reduced 
when increasing the saturation induced by the probe laser 

In order to expose the physical mechanism responsi- 
ble for the loss of coherence as clearly as possible, we 
will consider in this paper two two-level atoms in free 
space - the simplest system exhibiting coherent backscat- 
tering. Effects which arise in the presence of a larger 
number of atoms, such as the nonlinear index of refrac- 
tion of an atomic medium, will be relegated to future 
publications. Furthermore, we assume that the distance 
between the two atoms is much larger than the optical 
wavelength, such that we can neglect recurrent scattering 
(corresponding to a dilute medium in the case of many 
atoms). To calculate the photodetection signal of the 
light emitted by the two atoms, we use scattering the- 
ory. Generally, the higher the intensity of the incoming 
light, the more photons are scattered inelastically. In 
the present paper, we restrict ourselves to two-photon 
scattering. Thereby, we employ a perturbative approach, 
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valid up to second order in the incident intensity. 

The paper is organized as follows. In Sec. II, we sum- 
marize known results about the scattering of two photons 
by a single atom. After introducing the scattering opera- 
tor in Sec. II B, we obtain the corresponding photodetec- 
tion signal in Sec. II C, thereby recovering the resonance 
fluorescence spectrum in second order of the intensity. In 
Sec. Ill, we add a second atom to our model. We pro- 
ceed in a similar way as in Sec. II, using the results of the 
single-atom case as a building block of the two-atom solu- 
tion. After deriving the scattering operator in Sec. Ill A, 
we calculate the photodetection signal in Sec. IIIB. In 
contrast to the single-atom case, the latter contains in- 
terference between the light emitted by the two atoms, 
enhancing the detection signal in the backscattering di- 
rection. In this way, we obtain the main result of the 
present paper, the backscattering enhancement factor, 
which is found to be smaller than two, due to inelastic 
scattering. This fact is interpreted in Sec. Ill C as a loss 
of coherence between the light scattered by both atoms 
in opposite order. Regarding the undetected photon as a 
path detector for the detected photon, we can explain the 
loss of coherence by an analogy to the double-slit experi- 
ment, where the interference pattern is washed out if we 
try to observe which slit the particle has passed through. 
Finally, Sec. Ill D concludes the paper. 



II. SINGLE ATOM 

Let us start with discussing the scattering of two pho- 
tons by a single atom. This is useful since we will assume 
later that the second atom is far away from the first one. 
The two-atom scattering process can then be viewed at 
as a succession of two single-atom scattering processes. 



A. Approximations and Hamiltonian 

We assume a two-level atom located at a fixed position 
r. As already mentioned above, neglecting the exter- 
nal atomic motion is justified at very low temperatures, 
where the Doppler shift induced by the atomic motion is 
small enough. Also the recoil effect, i.e. the change of 
the atomic velocity when scattering a photon, can be ne- 
glected - provided that the number of scattering events is 
not too large. On the other hand, the temperature should 
still be high enough such that the external atomic mo- 
tion need not to be treated quantum mechanically. Fur- 
thermore, let us stress that we consider an undegenerate 
atomic ground state (J = 0). This is important since 
coherent backscattering may be severely affected by de- 
generacy [11-14]. The excited state is then threefold de- 
generate (J = 1). Which one of the three excited states 
is populated depends on the polarization of the absorbed 
photon. 



With the approximations mentioned above, our Hamil- 
tonian reads as follows: 



H = H Q + V, 

where 



(1) 



(2) 



k,s 



V = E (^^V • e ks )a ks - rge-^io ■ e ks )4 s ) , (3) 

k,s 

denote the free evolution, and the interaction, respec- 
tively (in units where % — 1). Here, the operators and 
a describe transitions between the atomic ground and ex- 
cited states, with energy difference cj at (in the case of an 
isolated atom), whereas aj^ and a^ s create and annihi- 
late a photon in mode k (a plane wave with wavevector k) 
and polarization eks (perpendicular to k). The coupling 
constant 



Wk 

2e L3 



1/2 



d 



W a t 

2e L 3 



1/2 



(4) 



with L 3 the quantization volume (which will finally drop 
out of the equations, when taking the limit L — ► oo) and 
d the magnitude of the atomic dipole, determines the 
strength of the atom-field coupling. 

In Eq. (3), we have employed the so-called 'rotating 
wave approximation': a transition from one of the excited 
states to the ground state is only possible by emitting a 
photon, and vice versa by absorption. This is justified 
since we will restrict ourselves to near-resonant processes, 
where only photons with frequencies close to the atomic 
resonance are important (i.e. |wk — w a t| <C w at ). For 
the same reason, we may assume a constant value of g in 
Eq. (4), i.e. neglect its dependence on uj^. 

Due to the coupling to the electromagnetic vacuum, 
the state |e) is unstable: after an average lifetime given 

by 
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(5) 



an excited atom decays into the ground state, through 
spontaneous emission of a photon. This gives rise to an 
effective, complex atomic resonance frequency 



= Wat 



'2' 



(6) 



where also the real part w at is shifted, as compared to 
the isolated atom, Eq. (2). 



B. Scattering matrix 

In the following, we make use of scattering theory in 
order to calculate the properties of the light emitted by 
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the atoms. Here, the object of interest is the scattering 
operator S, which connects the initial and final photon 
states |z) and |/): 

l/> = 5|i). (7) 

The initial and final state of the atom is always the 
ground state \g), which we do not explicitly write in the 
following. Furthermore, we will restrict ourselves to the 
scattering of two photons, thereby employing a pertur- 
bative approach, valid up to second order in the incident 
intensity. 

Since, as we will see below, the two photons may be 
scattered independently from each other, we consider 
first the scattering of a single photon. In order to dis- 
tinguish between the scattered and non-scattered part 
of the photon wavepacket, the transition operator Ti is 
introduced as follows: 



Si = 1 - 2m8(uj f -Wi)Ti, 



(8) 



where the 5-function implies conservation of the photon's 
frequency (which follows from energy conservation, since 
the state of the atom is the same before and after scat- 
tering). For one-photon states, its matrix elements read 
[17]: ' 

9 2 



(k/e/lTilk^) = 



■ w 



(d ■ e))e i < k '- k '>- r '. (9) 



The situation changes when considering a second pho- 
ton. It is convenient to write the matrix elements in the 
following form: 1 

(k 3 e 3 ,k 4 e 4 |,S' 2 |k :l ei,k 2 e 2 ) = 

(k 3 e 3 |S' 1 |k 1 e 1 )(k4e4|5 , 1 |k2e2) 

+ (k 3 e 3 1 | k 2 e 2 ) (k 4 e 4 1 Si | ki ei ) 
+<k 3 e3,k 4 e 4 |T 2 |k 1 ei,k 2 e 2 }. (10) 

Here, the first two terms scatter the two photons indepen- 
dently from each other. (There are two terms since the 
photons are indistinguishable: the final photon |k 3 e 3 ), 
for example, may correspond either to the initial photon 
|ki£i) or |k 2 e 2 ).) Since, however, the atom cannot inter- 
act with the second photon while it is excited by the first 
one, the photons are in fact not completely independent. 
This gives rise to the second term [18]: 

(k 3 e 3 ,k 4 e4|T2|k 1 e 1 ,k 2 e 2 ) = 
.g 4 S(u>i + u)2 — u>3 — w 4 ) 



2ni 



+ 



(u>l - U> ){u>2 - W ) VW3-W0 W 4 -W 0/ 

[(ei • e* 3 )(e 2 ■ e* 4 ) + (e 2 • e$)(ei • e* 4 )}e i ^ +k ^-^ r . (11) 



Although their sum is conserved, the individual frequen- 
cies of both photons may be changed by T 2 , for what 
reason we call it 'inelastic' scattering. 



C. Photodetection signal 

Given the final photon state |/), the intensity of the 
photodetection signal, as measured by a broadband de- 
tector (polarization cd) located at R at time t reads [17]: 



(/|25(->(R,t)25<+>(R,t)|/>. 



(12) 



Here, the detection of the photon is described by the 
electric field operator 



(R, t) = I ]T(6 ks • eBK^-Vs, (13) 

k,s 



which annihilates a photon at position R. 

As initial state, we consider a state of N photons 

|*jv> = V^V! h(ki)...ft(kj V )|ki€i,...,kj V e L ), 

(ki...kjv) 

(14) 

where all photons are described by the same single- 
photon wavepacket 



h)=J2h(k)\ke L ). 



(15) 



The factor y/N\ in Eq. (14) arises from the symmetry 
under exchange of photons as bosonic particles, and is 
required to obtain the correct normalization, (zat |*tv) = 
(*i = 1- We assume that the wavepacket describes an 
almost plane wave, i.e., h(k) is sharply peaked around 
its center k^ ('sharply' means: much narrower than T). 
For this reason, we may also neglect in Eq. (14) the de- 
pendence of the initial polarization vector el on k. 

The initial state |ijv) corresponds to the following in- 
cident intensity seen by the atom at position r and 
time t = 0, obtained by inserting |ijv) instead of |/) in 
Eq. (12), and summing over the detector polarization e D : 



L - - N ^f 



^> ik -Mk) 



(16) 



In the following, we will use a dimensionless quantity, the 
so-called 'saturation parameter' 



1 Eq. (10) is valid only if kiei / k 2 £2 and k 3 e 3 / k 4 e 4 . We 
will not consider double occupancy of modes in the following, 
since it can be neglected in the continuous limit of infinite 
mode density. In other words: two photons are never exactly 
in the same mode, although they may be infinitesimally close 
to each other. 



2d 2 L r 



2Ng* 



\U>L ~ W 



\u L - w 1 



5> ik - r Mk) 



(17) 



It accounts for the fact that photons interact less strongly 
with the atom if they are far detuned from the atomic 
resonance (i.e. if \u>l — ojq\ is large). From the solution of 
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the optical Bloch equations [17], it is known that s deter- 
mines the ratio between inelastic and elastic scattering, 
see Eq. (37) below. 

We are interested in the photodetection signal mea- 
sured at position R at time t = |R — r| (the time needed 
for the scattered light to reach the detector, in units 
where c = 1). We assume that the detector is placed far 
away from the atom, such that t = |R — r| is long enough 
for the scattering approach to be valid. Furthermore, the 
detector should not be placed in the direction of the ini- 
tial wavevector k^, such that only scattered photons are 
detected (i.e. P<+) (R, t)|ii) =0). 

In order to proceed, we have to generalize the scat- 
tering operator for two photons, Eq. (10), to the case of 
N photons. For this purpose, we assume that the satu- 
ration parameter s is so small that at most one photon 
pair is scattered inelastically. This yields the photodetec- 
tion signal up to second order in s, see below. Summing 
over the different pairs and taking into account all 

possible permutations of the N photons, we obtain: 

N 

(k' 1 ...k / JV |5 JV |k 1 ...k Ar )=^n^i s ii k ^w> + 



with 



JV-l 



(ki...k JV _i|V>i) = ViVLE IJ<ki|/i) 



(ki . . .kjv-ilV^) 



i=i 

N-l 



N-l 



£<kiisi> n< k *i/i>. 



(23) 



(24) 



»=i 

JV-l 



i=i 



N-l 



i,j = l 

i<j 



1=1 
l&,3 



and 



£ = <0|£ (+) (R,*)|/i>, 
\ gi )=E(+\R,t)\g 2 ). 



(26) 
(27) 



+ 



N 

E E 

i,j = l P N /P2 

i<j 



According to Eq. (12), the norm I = (ip\tp) gives the total 
intensity. Let us first concentrate on the contributions 
from and \i>2)- (As we will argue later, ^3) can be 
neglected). We obtain a sum of three terms, from elastic 
and inelastic scattering, and their interference. Using 
= 1 (since S\ is unitary), we obtain: 



(k^.|T 2 |kp w(<) k Pw(i) > n (k;i^|k PN(0 ). 



1=1 



(In the following, we do not write explicitly the polariza- 
tion vectors.) Eq. (18) contains a sum over all permu- 
tations Pjv of the N indices {1 . . . N}, modulo a permu- 
tation of the two indices Pjv(«) and Pjv(j) in the second 
term, where the latter permutation is included in the 
two-photon operator T 2 , see Eq. (11). In the case N = 2, 
the above expression agrees with the one of the previous 
section, Eq. (10). According to Eq. (18), the final photon 
state |/jv) = Sjvlijv) can be expressed as follows: 

N 



(18) = {^) = n\e\ 2 , 

= N(N-l)Re{V2E(g 1 \f 1 )}, 
I in = (V, 2 |V> 2 ) = % i\9i) 

+ n{n -T- 2) \«m\>. 



(28) 
(29) 

(30) 



(k 1 ...kj V |.fjv> = VAnn< k 'l/i> + 

1=1 

Aah n n 

+ y— e^i^> n <kziy!>, 



(19) 



i,j=l 

i<j 



1=1 



in terms of the one- and two-photon states 

l/i>=Si|ii>, 
|52> =T 2 \i2). 



(20) 
(21) 



Following Eq. (12), we now apply the electric field op- 
erator on the final photon state. It may annihilate either 
an elastically or an inelastically scattered photon. Corre- 
spondingly, we obtain the following three contributions: 



fi 2) 

Whereas in P^ ' , the frequency of the detected photon 
is fixed to u> — ujl (since one-photon scattering is elas- 
tic), this is not the case for L m , where the overlap (<7i|<?i) 
implies an integral over cu. Thereby, we obtain an elastic 
and inelastic component of the detection signal. 

To complete the calculation, we insert the one- and 
two-photon scattering matrices given in Sec. II B. Using 
Eqs. (8,9), the final one-photon state reads: 

l/i) - *i - — x 

LUl — LJq 

x Yl /i(ki)^ i -o; / )(e i -e}) e ^- k /)' r |k / e / ). (31) 

ki,k f e f 

Since the wavepacket h(ki) is quasi- monochromatic, we 
may replace the argument of functions which vary slowly 
(i.e. on the scale of T) as a function of by the constant 
value uil- Applying the electric field operator on see 
Eq. (26), yields, under the assumptions given above: 



= 25(+>(R,t)|/>= 



(22) 



E - -3r(e L ■ e* D )g ^ J^uiy.. 



(32) 



i=i 
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Similarly, we obtain for the inelastic part, see 
Eqs. (11,21,27): 

2 ^ V2^ k -- k /)^(k 4 )(£L-e}) |1 . 



(w/ - w )(2w L - w/ - ujq) 



(33) 



According to Eqs. (28-30), we obtain the following inten- 
sity I = 1^' + 1^' + I[ n of the photodetection signal: 

^=4 4 2) = V^ 2 , (34) 



N-l 



N 2 



N 



with the prefactor: 



V 



f3T\e L 



V AdLj L R 



(35) 



(36) 



The term proportional to |(/i|gi)| 2 in Eq. (30) gives a 
contribution to the inelastic component in third order of 
s, which can be neglected. As it should be, for large 
N - such that the first photon can be absorbed without 
significantly changing the saturation induced by the re- 
maining (N — 1) photons - the above result agrees with 
the elastic and inelastic components 



Id = V 



2(1 + s) 2 ' 



Iin = V 



2(1 + .s) 2 



(37) 



of the resonance fluorescence as predicted by the Bloch 
equations [17], expanded up to second order in s. 

However, we have not yet accounted for the third term 
\fa) in Eq. (22). If we compare Eqs. (19,23,25), we note 
that \fa) + {ips} = >/~NE\fN-i), and hence the norm of 
\fa) + \fa) equals the norm of \fa), Eq. (28), provided 
that the norm of |/iv-i) is 1. Although the latter condi- 
tion is not necessarily fulfilled if the scattering operator 
is truncated as in Eq. (18), its unitarity will be recovered 
when including higher scattering orders. Similarly, it can 
be shown that contributions from (fa \fa) - if they are not 
of third order in s - are exactly canceled by other terms 
which appear in (fa\fa) when including into (fa\ another 
inelastically scattered photon pair. Hence, the term ^3) 
does not contribute to the photodetection signal up to 
second order in the saturation parameter s. 

By putting a spectral filter in front of the detector, 
we can resolve the power spectrum P(u>) of the detection 
signal, i.e. the probability of detecting a photon of a 
definite frequency uj. Since elastic scattering conserves 
the frequency, the spectrum exhibits a sharp peak at wi 
(almost a (5-function for our quasi-monochromatic initial 
wavepacket /), 



P(w) =I {cl) 6 f {u-Lu L )+p( in \u), 



(38) 



whereas the inelastic component depends smoothly on uj. 
The latter is proportional to the absolute square of the 



P ,m) (w) 

(i (m) /r) 
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FIG. 1. Inelastic resonance fluorescence spectrum p( m ' (a;), 
Eq. (40), for small saturation, s <C 1, (a) zero detuning 
S — u>l — Wat = 0, and (b) S = 2Y. The dashed lines in- 
dicate the position of the elastic peak at lul, see Eq. (38). 

inelastic transition amplitude, Eq. (11) (with wi = 
u>2 = <jJl the initial frequency, u>3 — u> the frequency 
of the detected photon, and uu^ = 2ul)l — u). With the 
correct normalization, 



we obtain: 



p (in) H = 



r/( in ) 

47T 



1 



+ 



u> — Wo 2u>l — uj — luq 



(39) 



(40) 



For zero detuning, uj l — w at , the inelastic spectrum con- 
sists of a peak of width 0.64r, whereas for large detuning 
S = ujl — uj at (he. if 4<5 2 3> T 2 ), there are two peaks of 
width r at uj = ujl ± 8, see Fig. 1. 2 Note that one of 
them is centered exactly at the atomic resonance. Evi- 
dently, this will be important if we allow the scattered 
photons to interact with a second atom, as we will do 
now. 



III. TWO ATOMS 

A. Scattering matrix 

Let us now turn to the case of two atoms alone in 
vacuum. We assume that the second atom is far away 
from the first one, compared to the optical wavelength. 
This means that we may restrict ourselves to processes 
where at most one of the two photons is scattered by both 



2 The reader may have in mind that the resonance fluores- 
cence actually exhibits three peaks [19]. However, the one at 
u) — u)l is of higher order in s, since it arises from three- 
photon scattering. 
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atoms. As shown in appendix A, the corresponding scat- 
tering matrix can then be obtained in a simple way from 
the single-atom scattering matrix, see Eqs. (A9,A10): 
apart from the geometrical phase factors e ±4k ' ri < 2 for ab- 
sorption or emission of a photon |k) by atom 1 or 2, and 
the terms depending on the polarization, we only have to 
take into account the 'photon exchange factor' 



B{u>) 



ZTe 1 



Aujr 12 (uJ - wq) ' 



(41) 



depending on the frequency u) of the doubly scattered 
photon. In Eq. (41), we recognize the propagation of a 
spherical wave from one atom to the other one (inversely 
proportional to their distance ri 2 ), and the amplitude 
(w — Ldo)^ 1 describing scattering by a single atom, see 
Eq. (9). 

In particular, the structure of the scattering opera- 
tor as a sum of an elastic single-photon and an inelastic 
two-photon component is the same as before, compare 
Eqs. (8,10): 

(2) 

(k 3 e 3 , k 4 e 4 |S^ ; |kiei,k 2 e 2 ) = 

(k 3 e3|5p ) |k 1 e 1 )(k 4 e 4 |5f ) |k 2 e 2 ) 

+ (k 3 e 3 |S 1 2) |k 2 e 2 )(k 4 e 4 |Sf ) |k 1 e 1 ) 
+ (k 3 e 3 ,k 4 e4|T 2 (2) |k 1 e 1 ,k 2 e 2 >, (42) 

where the single-photon component of S 1 - 2 ^ contains also 
the non-scattered wave, see diagram (la) in Fig. 2: 



S[ 2) = 1 - 2m5(uj f ~uo i )T[ 2) . 



(43) 



The remaining single-photon processes are also shown in 
Fig. 2. The photon may be scattered by only one atom (1 
or 2), or by both (first 1, then 2, and vice versa). Corre- 
spondingly, the single-photon transition operator reads: 

( k/e/ |Tf |k iei ) = -^-x( 

gi^-k,).^ ^ . e * } + B ( w .)( e . . Al2 . c *) e <M**-'i>) 

e i(k 4 -k /) .r a ( (e . . e * } + B ( w .) (e . . Al2 • e })e iki '( r2 - ri )) | 

As mentioned above, for the two double-scattering pro- 
cesses, see diagram (Id) and (le) in Fig. 2, we have to 
multiply the one-atom transition operator (k/e/|T|kje,), 
Eq. (9), with the photon exchange factor B(u>i), see 
Eq. (41), and to adjust the geometrical phase factor. Fur- 
thermore, the fact that the photon propagates in the di- 
rection r 2 — ri between the two scattering events implies 
a projection Ai 2 of the polarization vector onto the plane 
perpendicular to r 2 — ri. Thereby, the term a ■ (for 
scattering by a single atom) is replaced by ej • Ai 2 • e^. 

In the case of inelastic two-photon scattering, the dou- 
bly scattered photon may be scattered first inelastically 



(by atom 1 or 2), and then clastically (by the other atom), 
or vice versa, compare, e.g., the diagrams (2a) and (2d) 
in Fig. 3. Correspondingly, the frequency to be inserted 
in the photon exchange factor B(u), Eq. (41), is either 
the final or initial frequency of this photon, see Eq. (A9) 
or Eq. (A10). In total, we obtain: 



(k 3 e 3 ,k 4 e4|T 2 (2) |k iei ,k 2 e 2 ) = 



2iri 



.g 5{oj\ + lu2 



1 1 

+ 



(Wl - W )(w 2 - O> ) \W3-UJ0 LUi-LUo 

e i(k 1+ k 2 -k3-k 4 ). ri x |( ei . e *)( £2 . e *) 

+ B(wi) (ei ■ A12 ■ e3)(e 2 ■ ej) e^-^-rO 
+ B{lo 2 ) (e, ■ e* 3 )(e 2 ■ A 12 ■ e* 4 ) e *^-*0 
+ B(uj 3 ) (ei • A12 • e 3 )(e 2 ■ el) e^'^-n) 

+ B(ui) (ei • e*)(e 2 • A 12 • e%) e -*4-(^-'0 

+ e i(ki+k2-k 3 -k 4 ).r 2 x I (£i . £ * )(£2 . £ *) 

+ B(ui) (ei • Ai 2 • el){e 2 ■ e\) e*-^"^) 
+ B(co 2 ) (e 1 -e*)(e 2 -A 12 . e *) ^^-^) 
+ B(uj 3 ) (ei • A12 • e 3 )(e 2 • el) e"^'^-^) 



+ B(ui) (ei • e*)(e 2 • A 12 • e* 4 ) e"*-^-^) 



+ ^kiei <-> k 2 e 2 ^) 



(45) 



Here, the last line denotes additional terms arising from 
exchanging the initial (or, equivalently, final) photons. 
We recognize two terms describing the scattering by atom 
1 or 2 alone, see diagrams (2i) and (2j) in Fig. 4, and 
eight different terms describing the processes where both 
atoms are involved, see diagrams (2a-2h) in Figs. 3 and 
4. Note that the terms depending on the polarization 
allow to identify the photon which is scattered by both 
atoms. This photon is marked by full arrows in Figs. 2-4, 
whereas the open arrows denote the photon scattered by 
(44)only one atom. If we assume that |k 4 e4) is the doubly 
scattered photon, the ten terms in Eq. (45) correspond 
(from top to bottom) to the diagrams (2i), (2e), (2a), 
(2h), (2d), (2j), (2f), (2b), (2g), and (2c), respectively. 



B. Direct calculation of the enhancement factor 

Having at hand the scattering matrix, we now deter- 
mine the intensity of the photodetection signal. In prin- 
ciple, the calculation can be performed in the same way 
as in the single-atom case, Sec. II C. However, the detec- 
tion signal will now depend non-trivially on the position 
of the atoms and the detector, due to the fact that 
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(Id) (le) 




FIG. 2. Scattering of a single photon by two distant atoms. 
In the coherent backscattering experiment, only the doubly 
scattered photon is detected, see diagrams (ld,e). Conse- 
quently, the diagrams (la-c), with single or no scattering, 
describe the undetected photon. 




FIG. 3. Inelastic scattering of two photons by two distant 
atoms. Only the doubly scattered photon (full arrows) is 
detected. Since the photon frequencies are changed by the 
inelastic scattering event at the atom where both photons 
meet, the amplitude of the elastic scattering event at the sec- 
ond atom depends on whether the inelastic scattering occurs 
before the elastic one (2c, d) or after (2a,b). 

the photons emitted by one atom interfere with the 
photons emitted by the other one. Even if we average 
over the positions ri and r2 of the atoms, the interference 
is not completely washed out. There remains an 




FIG. 4. Remaining diagrams describing inelastic scattering 
of two photons by two atoms. In the coherent backscattering 
experiment, they are filtered out by using the h \\ h polar- 
ization channel (see Sec. IIIB), in which a singly scattered 
photon (open arrows) cannot be detected. 

enhanced probability to detect a photon in the direc- 
tion opposite to the incident wave - an effect which is 
known as coherent backscattering. In the case of two 
atoms, it arises from double scattering: in the backscat- 
tering direction, a photon scattered first by atom 1, and 
then by atom 2, intcrfers constructively with the corre- 
sponding reversed path. 

In order to examine cleanly this interference effect, we 
therefore assume that only doubly scattered photons are 
detected. Experimentally, this can be realized by using 
circularly polarized light = (l,i,0)/v2 (in Euclidean 
coordinates, where the z axis is parallel to k^) , and de- 
tecting the scattered photons in the helicity preserving 
channel = e* L . This implies (ej, • e* D ) = 0, i.e. no 
singly scattered photons can be detected in the helicity 
preserving polarization channel. If we look at the in- 
elastic part of the scattering matrix, Eq. (45), assuming 
(without loss of generality) that the photon |k4e4) is the 
detected one, this means that all terms with (ei • e|) or 
(e2 • e|) are filtered out. These are the diagrams shown 
in Fig. 4, and only those of Fig. 3 remain. Concerning 
the elastic single-photon scattering, see Fig. 2 we keep 
the single-scattering diagrams (la-c) to describe the un- 
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detected photon. For the sake of completeness, we will 
repeat in appendix B the following calculation for the 
case of scalar photons, where, a priori, all the diagrams 
shown in Figs. 2-4 contribute. 



1. Elastic contribution 

Let us begin with the contribution 1^ of one-photon 
scattering. According to Eqs. (26,28), it is obtained by 
applying the electric field on the final state of single- 
photon scattering. As explained above, only the dia- 
grams (Id) and (le) in Fig. 2 contribute. At first, we 
concentrate on the phase factors depending on the posi- 
tion of the atoms. If is the wavevector of the incident 
photon, and the detector is located in the direction kjj 
(with \ko\ = |ki|, since one-photon scattering conserves 
the frequency), we obtain exp(iri • k^ — ir2 ■ k^) for (Id) 
and exp(ir2 • k^ — ir\ • kjj) for (le). Evidently, the phases 
are identical if kp = — k^, i.e., (Id) and (le) interfere 
constructively in the backscattering direction. On the 
other hand, if k^ ^ kp (more precisely: if the angle be- 
tween k^ and is much larger than some characteristic 
quantity 8c), the interference between (Id) and (le) van- 
ishes when averaging over the positions of the atoms. For 
simplicity, we fix the distance r\2 and average only over 
the angular variables of ri — Y2- In this case, the width of 
the enhanced backscattering signal (which is also called 
'the cone') is given by 9q — 1/ (wr 12 ). In total, we obtain 
both for the background intensity (known in the litera- 
ture as the 'ladder term' L), and the additional intensity 
in backscattering direction (the 'crossed term' C) twice 
the result r]s /2 of the single-atom case, 



L (D =C (i) 



T]S, 



apart from a modification of the prefactor 

2 

/ .51 \ 

V 

2 



(4ib) 2(|BK)|2|e - Al2 ' e - |2> - 



9T Z 



lQ<L) L Rr\2\uL — wo | 



(46) 

(47) 
(48) 



Here, Eq. (47) implies an average over the positions of 
the two atoms. The polarization-dependent term \er, ■ 



c * |2 



sin 6/4 is given by the angle 9 between the 



incident laser k^ and the two atoms V\2 = T\ — r 2 . Then, 
a spherical distribution of Y\2, at fixed distance r\2, yields 
the result given in Eq. (48). The fact that = C^> 
can be traced back to the reciprocity symmetry [4]. 

Next, we examine the interference between two-photon 
and one-photon scattering, which gives rise to the elastic 

(2) 

component of the intensity in second order of s, see 
Sec. IIC. According to Eq. (29), 1^ is given by the 
overlap of the respective quantum states |/i) and \g\) 
of the undetected photon, which amounts to a sum over 
the latter's state |ke) (i.e., (gi\fi) = E kie (#i|ke)(ke|/i)). 



First, we concentrate on the phase factor exp(— ik • ri^) 
of the undetected photon, depending on whether it is 
emitted by atom 1 or 2. Integrating over the angular 
variables fl^ of k (at fixed |k| = ojl), we obtain, if |k) is 
emitted by different atoms: 



/ 



rfO k e^-3) = 47 r Sin(a;Lri2) « 1. (49) 

ULT\2 



Since we have assumed wl^12 S> 1, the above term can 
be neglected. In other words, diagrams where the unde- 
tected photon is emitted by different atoms do not inter- 
fere in leading order of l/(w L r 12 ). If we now select one of 
the four diagrams (2a-d) describing two-photon scatter- 
ing, we can discard among the three one-photon diagrams 
(la-c), the one where the undetected photon is scattered 
by the 'wrong' atom. The remaining two exactly give the 
final state of a photon scattered by a single atom, as de- 
scribed by Eq. (20). Concerning the detected photon of 
the one-photon scattering, we can choose either diagram 
(Id) or (le). As already discussed above, one of them 
gives a contribution to the background L, and the other 
one to the enhanced backscattering signal C. As there 
are in total four diagrams (2a-d), we obtain both for L 
and C four times the result —rjs 2 of the single-atom case: 



L (2,el) = c (2,cl) 



-4?7s 2 , 



(50) 



Note that the total elastic ladder term, Eq. (46) and 
Eq. (50), equals the total elastic crossed one, Eq. (46) and 
Eq. (50). This means interference with maximal contrast, 
corresponding to the maximal possible enhancement fac- 
tor of two. 



2. Inelastic contribution 

The inelastic component /; n of the intensity, finally, 
arises from two-photon scattering. Here, the overlap 
(<?i|<?i), see Eq. 30, again implies a sum over the un- 
detected photon, which now may have a frequency dif- 
ferent from ujl- With two atoms, \gi) is a sum of 
four different contributions, corresponding to the dia- 
grams (2a-d) . Correspondingly, we obtain diagonal terms 
(|2a| 2 , . . . , |2d| 2 ), which contribute to the background sig- 
nal, and interference terms, which may contribute to the 
backscattering cone, see below. 

Let us examine first the diagrams (2a) and (2b), where 
the elastic scattering event occurs before the inelastic 
one. Here, the single-atom scattering amplitude is mul- 
tiplied by a constant factor B(lul). This means that - 
apart from the modification of the prefactor rj - both 
|2a| 2 and |2b| 2 give the same result as in the single-atom 
case, Eqs. (30,35): 

I II = lJdu,pW{u,) = ij£ (51) 
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In the other two cases (2c) and (2d), the frequency to be 
inserted in the factor B(uj), Eq. (41), equals the final fre- 
quency of the detected photon (or - equivalently - of the 
undetected one, since |-B(w)| 2 = \B(2lul — u)\ 2 ). Hence, 
a factor |£(u;)| 2 must be inserted in the integral over the 
inelastic power spectrum, Eq. (39). The resulting inte- 
gral can be easily performed, and yields: 



= 1 
V 

„s 2 



du> 



LU L - UJq 



LU — LUq 

S 2 



3 

4 + f2 



(52) 



Hence, the four diagonal terms |2a| 2 , . . . , |2d| 2 , give the 
following contribution to the inelastic background inten- 
sity: 



L (in) = 2Ji + 2J„ = 



(53) 



Note that, for 8 — 0, the contribution from (52) is smaller 
than the one from (51) (by a factor 3/4). This is due to 
the fact that, after the inelastic scattering event, the pho- 
ton frequencies are no longer exactly on resonance, see 
Fig. 1(a), what reduces the cross section of the scattering 
by the other atom. The opposite is the case for large de- 
tuning 8: here, the inelastic scattering brings one of the 
two photons close to the atomic resonance, see Fig. 1(b), 
thereby increasing the corresponding contribution to the 
background signal. 

The inelastic component of the enhanced backscatter- 
ing signal arises from the interference of (2a) with (2d), 
and (2b) with (2c). (Remember that every diagram in- 
terferes only with those where the undetected photon is 
emitted by the same atom.) As argued above, equality 
of the corresponding geometrical phases, and thereby full 
constructive interference, is guaranteed if the wavevec- 
tor of the detected photon is opposite to the incident 
wavevector, i.e., ko = — k^. Obviously, this condition 
will not be exactly fulfilled in the presence of inelastic 
scattering, even in exact backscattering direction (since 
in general |ko| ^ |ki|). The difference can be neglected, 
however, if we assume that the atomic linewidth T and 
the detuning 5 = uol — w a t, i.e. the parameters which de- 
termine the width of the power spectrum, see Fig. 1, are 
much smaller than the inverse of the distance r\i between 
the atoms: 



<5,T< 



ri2 



(54) 



In other words: the propagation time ri 2 /c between the 
atoms is much smaller than the timescales associated 
with 8 and T. This condition ensures a vanishing geo- 
metric phase difference, i.e. exp[(ki + k£i)-(ri— r 2 )] — 1, 
and is well fulfilled in the experiment [16]. What remains 
is the integration over the inelastic spectrum, taking into 
account the photon exchange factors B(lol) or B(u>) in 
the cases (2a, b) or (2c, d), respectively: 



2 [<L;Re{^—^)pW(u 
J ( w-wq J 



-rjs 2 . 



(55) 



Here, we have neglected the exponential factor 
e 2(w-w L )ri2 r*j i describing the propagation in the vac- 
uum - the same approximation as above, see Eq. (54). 
From the two interfering pairs of diagrams, the inelastic 
contribution to the backscattering signal is obtained as 
twice the result of Eq. (55) with modified prefactor fj\ 



C (in) 



(56) 



Note that is strictly smaller than the inelastic back- 
ground, Eq. (53), which leads to a reduction of the 
backscattering enhancement factor, see below. This is 
consistent with the fact that two interfering diagrams, 
e.g. (2a) and (2d), are no more linked by the reciprocity 
symmetry: only diagrams with identical initial and fi- 
nal photon frequencies interfere with each other, whereas 
the reciprocity symmetry connects diagrams where initial 
and final frequencies are exchanged. 



3. Double scattering enhancement factor 

Adding all contribution, we have: 

L = L« + + =ij^8- Is 2 + ^ (57) 

C = + + =rj(s- |s 2 ) . (58) 

Finally, the double scattering enhancement factor reads: 
L + C 8 - (19 - 4S 2 /T 2 )s 



a = 



4- (9-45 2 /r 2 )i 



(59) 



Remember that single scattering has been removed by 
the helicity-preserving polarization channel. 

At this stage, it is convenient to introduce the satura- 
tion parameter on resonance: 



so 



2d 2 I in 

r 2 /4 



AS 2 " 

1 + f2 I 



(60) 



which depends only on the incident intensity I ln (and not 
on the detuning 6). Then, Eq. (59) can be rewritten: 



with 



2 + .T 

l + x' 



so s 



4 -10s 4 



(61) 



(62) 



Here, we have used that s is small - otherwise, our per- 
turbative treatment (two-photon scattering) would be in- 
valid. If the detuning 8 is of the order of the linewidth 
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FIG. 5. The enhancement factor a — (8 + so)/(4 + so) as 
a function of the incident intensity so, and large detuning 
8 — lol — w a t- If 8 is not large, the displayed curve is valid 
only up to intensities so <C 1 + 4<5 2 /T 2 , corresponding to a 
small saturation parameter s <C 1, cf. Eq. (60). 

T, this implies that so must also be small. In this case, 
Eq. (61) yields: 



so 



(63) 



In principle, however, we may choose also a large value 
of the detuning 6, as long as we stay near-resonant, and 



fulfill (5/r) 2 < l/(wri2 



This means that sq may be 



large although s is small, see Eq. (60). In that case, the 
enhancement factor is given by Eq. (61), with x = so/4, 
see Fig. 5. This equation is valid for all values of s 
corresponding to small s, i.e. sq <C 1 + 4S 2 /T 2 . 

It may appear surprising that the enhancement fac- 
tor a depends only on the intensity so of the incident 
light, see Eqs. (61,62), whereas the intensity scattered by 
a single atom is determined by the saturation parame- 
ter s, see Sec. II C. This result is related to the form of 
the inelastic spectrum, see Fig. 1: since one of the two 
photons is always close to the atomic resonance after the 
inelastic scattering, the asymmetry between the reversed 
paths (see the following section) is larger for larger ini- 
tial detuning 5, at a given value of s. Thereby, we can 
understand why, at a fixed s, the enhancement factor a 
decreases when increasing 5. However, we are not aware 
of an intuitive explanation why the relevant parameter 
turns out to be so, and not some other, similar combina- 
tion of 6 and s. 



3 This condition implies s \B(u a t)\ 2 < 1, see Eqs. (41,60), 
and thereby suppresses exchange of more than one resonant 
photon between the two atoms, leading to terms proportional 
to Sq (or higher order). 
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FIG. 6. Two reversed scattering paths, whose interference 
gives rise to enhanced backscattering. The left atom scatters 
inelastically, changing the photon frequency from ujl to ujd- 
Consequently, the amplitudes of the elastic scattering event 
by the right atom are different for both paths, see Eqs. (64,65). 
This leads to a reduced contrast of the interference pattern 
as a function of the detection angle 0, see Fig. 7. 



C. Interpretation: which-path information and 
coherence loss 

In this section, we discuss the physical mechanism re- 
sponsible for the reduction of the backscattering enhance- 
ment factor. As we have seen above, it originates solely 
from inelastic scattering. For this reason, we will only 
consider inelastic scattering in the following. 



1. Interference between counterpropagating paths 

Generally, coherent backscattering arises from con- 
structive interference between two scattering paths where 
the detected photon interacts with the respective scatter- 
crs in opposite order. The maximum enhancement factor 
of two is obtained if every path has a counterpropagating 
counterpart with the same amplitude. In the case of two 
photons, a 'scattering path' in principle also specifies the 
final state |k) of the undetected photon. As we have seen 
above, Eq. (49), the average over the angular variables 
of the undetected photon destroys interference between 
paths where the inelastic scattering occurs at different 
atoms (if the atoms are far away from each other) . Con- 
sequently, if we concentrate on the detected photon, we 
should compare only the two reversed paths where the 
inelastic scattering occurs at the same atom, and the fi- 
nal frequency lud = 2wl — is the same (due to energy 
conservation), as shown in Fig. 6. Here, the left atom is 
marked as the one which scatters inelastically. Neglect- 
ing the propagation in the vacuum, see Eq. (54), the am- 
plitudes Eiji of the two reversed paths are obtained by 
multiplying the scattering amplitudes of the elastic and 
inelastic scattering event, Eqs. (9,11). Since the elastic 
scattering occurs at two different frequencies u>l and lod, 
the amplitudes arc not identical: 
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Ei = 
En = 



1 



1 



U)D - t^O 
1 



+ 



2ui L - U) D 

1 



-ikr±0/2 



Wo / U>D — ^0 
ikr ± 6/2 



(64) I = j du D \Ei{uj D ) + En(u D )\' 



lu d - luq 2uj l -w D -u>oJ u L - uj 



(65) 



where 9 denotes the angle between detector and backscat- 
tcring direction, r± the perpendicular distance between 
the atoms, and prefactors not depending on or 6 are 
neglected. 

The asymmetry of amplitudes leads to a loss of con- 
trast in the interference pattern: 

\Ei + En\ 2 = \Ei\ 2 + \En\ 2 + 2\Ei\\Eii\cos{^ + kr^e), 

(66) 

if we plot the intensity as a function of the detection angle 
6, see Fig. 7(a). Here, the value of the phase shift reads: 



tan(</>o) 



2(6- e)T 
4<5e + r 2 ' 



(67) 



and — 7r < <f>Q < 7r, where e = lud — w a t is the detuning 
of the detected photon. The maximum contrast, where 
the intensity oscillates between zero and two times the 
mean value as a function of 6, is achieved only in the 
case = \Eu\. Nevertheless, also for \Ei\ ^ |£ii|, the 
two amplitudes interfere perfectly coherently: the con- 
trast is as large as it can be, given that the probabilities 
I-E-i.iiI 2 °f path I and II are different. This case is analo- 
gous to a double-slit experiment performed with a perfect 
monochromatic plane wave, but with different slit sizes. 
Here also, the contrast is reduced although the wave co- 
herence is perfectly preserved [21]. 



2. Second photon as a which-path detector 

What we have said so far is valid if the final frequency 
ujry is fixed. In reality, however, up is a random vari- 
able. This implies that the ratio Ei/Eu between the 
amplitudes of both paths fluctuates randomly, leading to 
a loss of coherence between the two paths. 4 If we look 
at the interference pattern of the average intensity, inte- 
grated over ojd, the loss of coherence reveals itself as a 
reduced contrast J max — / m ; n < A\fliln, i.e. smaller than 
the maximum value for two perfectly coherent waves of 
intensities h and In, respectively, see Eq. (66). The fac- 
tor 7i ji by which the contrast is reduced is called degree 
of coherence (see [21], p. 499-503): 



4 In general, fluctuations of the phase and of the absolute 
value of Ei /Eu both lead to 71,11 < 1, as defined in Eq. (69). 
In our case, the phase fluctuations have a stronger impact, at 
least for moderate values of the detuning 8 (not much larger 
than T). 



h + In + 2y/lihi 71,11 cos (4> + kr ± 6) , 



(68) 
(69) 



with Ii and In the average intensities of path I and II, 
respectively, given by Eqs. (51,52), and 4> the remaining 
phase shift. An example is shown in Fig. 7(b). To make 
the loss of coherence visible, we have indicated by the 
arrows the contrast obtained for the same values of Ji,n, 
but 71,11 = 1. 

An alternative physical explanation of the coherence 
loss can be obtained by an analogy to Young's famous 
double-slit experiment. As it is well known, interference 
is necessarily destroyed whenever we observe which slit 
the particle passes through (see, e.g., [22]). If \D{) and 
\Dn) denote the quantum states of the which-path detec- 
tor corresponding to path I and II, the degree of coher- 
ence is obtained as the overlap of the normalized detector 
states [20]: 



7iJi = 



KA|Ai)| 



^(DiIDiXAiIAi)' 



(70) 



This implies perfect coherence, 71,11 = 1, if the paths are 
indistinguishable (i.e. if the detector states are identical), 
and total loss of coherence, 71,11 = 0, if the paths can be 
distinguished with certainty (i.e. if the detector states 
are orthogonal). 

In our case, the path detector is given by the unde- 
tected photon. Remember that its frequency is correlated 
to the one of the detected photon, due to conservation 
of energy at the inelastic scattering event. Therefore, 
the different dependence of the amplitudes Ejjj of path 
I and II on the frequency of the detected photon, sec 
Eqs. (64,65), reflects itself in the final state of the unde- 
tected photon: 



\D I ) = J2E I (2LO L -w k )(e L -e*)\ke), (71) 

ke 

|Ai) = ^£n(2wL-w k )(e L -e*)|ke). (72) 

k£ 

Since \Di) and \Dn) are not identical, the state of the un- 
detected photon contains information about which path 
the first photon has taken. This leads to a loss of co- 
herence according to Eq. (70). Using Eqs. (64,65), wc 
obtain: 



7i,ii = 



9 + 4<5 2 /r 2 



12 + 16<5 2 /r 2 



tan(</>) = — , 



(73) 
(74) 



and — 7r < <j) < 7r. The probabilities of the two paths, 
h = (Di\D\) and In = {Dn\Du), respectively, have been 
calculated in the previous section, see Eqs. (51,52). Note 
that the phase shift, Eq. (74), vanishes for zero detuning 
- this can be traced back to the symmetry of the power 
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spectrum with respect to u) at , see Fig. 1(a), and the fact 
that the scattering amplitude (ujt, — Wo) _1 = — 2i/r is 
purely imaginary at resonance, u>l — ^at- 

One may wonder how the which-path information can 
be extracted from the second photon. In principle, this is 
possible, e.g., by performing a projection |£>i)(Di|, with 
possible measurement results '0' and '1'. If we measure 
'0', we know with certainty that path II has been taken, 
whereas in the case '1', both paths are possible, but with 
increased probability of path I. Due to the rather com- 
plicated expressions of \D{) and |-D n ), see Eqs. (71,72), 
it seems rather difficult to build such a detector in prac- 
tice. The simplest observation is to look at the frequency 
Wk = 2ojl — ojd- Since |^i(w_d)| 2 and \Eu(lud)\ 2 depend 
differently oiiwd, this observation does give some infor- 
mation about the path, although it does not resolve the 
phase dependence of E\{lod) and E\\{lod)- Nevertheless, 
the state of the second photon is completely determined 
by ojk (since the angular distribution is given by the po- 
larization, independent of I or II) . This implies that the 
which-path information can be erased by putting a spec- 
tral filter in front of the detector. Then, the measurement 
of the frequency gives no information (since it is always 
the same, determined by the filter), and, consequently, 
the coherence is fully restored, see Fig. 7(a). 



3. Coherence between the light emitted by atom 1 and 2 

At first sight, it may seem surprising that, for S > 0, 
the maximum of the interference pattern is not found in 
the backscattering direction, at 9 = 0. The reason is 
that we have specified the inelastically scattering atom, 
see Fig. 6, thereby introducing an asymmetry. However, 
with equal probability, the inelastic scattering may occur 
at the other atom, and then the interference pattern is 
shifted in the opposite direction. As already mentioned 
above, these two cases do not interfere, since they are 
distinguished by the undetected photon (in a similar way 
as the coherence between path I and II is reduced). If we 
add therefore both interference patterns incoherently, the 
new maximum is found at 9 — (as it should be), and 
the contrast of the interference pattern will be further 
reduced, by a factor cos(</>): 



2/i + 2/„ + 4V7iin 7i,ncos(0) cos(kr ± 9). (75) 



£ (in) 



£i(in) 



In the backscattering direction, 9 = 0, we recover the to- 
tal intensity as a sum of a background and an enhanced 
backscattering term and C^ n \ which we have cal- 
culated in the previous section, see Eqs. (53,56). 

What degree of coherence do we associate to the to- 
tal interference pattern? Note that we are now dealing 
with four different diagrams, see Fig. 3, whereas (second- 
order) coherence is a property of two interfering waves. 
Therefore, before we speak of coherence, we have to spec- 
ify how to divide the four diagrams into two waves. For 



that purpose, we consider the light emitted by atom 1 as 
the first wave, see diagrams (2a) and (2c), and the light 
by atom 2 as the second one, see diagrams (2b) and (2d). 
Indeed let us point out, even if it is obvious, that the total 
scattered electrical field is the sum of the electrical fields 
radiated by each atom in response to the local electrical 
field. This local field embodies the incident field and the 
field radiated by the other atom. This means that mul- 
tiple scattering is taken into account. Thus, the total ra- 
diated intensity contains interference terms between the 
fields radiated by each atom, some of them surviving the 
spatial and spectral averaging in the backscattering di- 
rection. In this respect, the CBS signal probes the spatial 
coherence between the two radiated fields. 

For reasons of symmetry, the intensities emitted by 
atom 1 and 2 are identical: 



h 



(76) 



According to the definition of the degree of coherence, 
hot =h+I 2 + ^Vhh 7i,2 cos(kr ± 9), (77) 

see Eq. (69), and taking into account Eqs. (75,76), we 
obtain 



71,2 = 7i.ii cos{(f>) . 

h + Jii 



(78) 



On the other hand, if we compare (78) with the under- 
braced terms in Eq. (75), the close relation between 71,2 
and the backscattering enhancement factor (consid- 
ering only inelastic scattering) becomes evident: 



Win) 

71,2 = T777T = a ' - f 



6 



7 + 4<5 2 /r 2 ' 



(79) 
(80) 



Here, we have inserted the results (53,56) of the previous 
section. From Eq. (78), we see that the coherence 71,2 
between atom 1 and 2 is reduced by the average over 
the power spectrum, leading to 71. n < 1, on the one 
hand, and by the random choice of atom 1 or 2 as the 
inelastically scattering atom, on the other one. The latter 
affects the coherence both due to different phases ±0 
and different amplitudes I\ and For not too large 
values of the detuning 5, the amplitude term can almost 
be neglected (being larger than 96% for S < T), whereas 
the phase term cos(</>), although strictly vanishing at S = 
0, gives a significant contribution if 5 is of the order of T, 
see Eq. (74). 

Let us note that the above interpretation in terms of 
a which-path experiment remains valid when adding the 
new pair of diagrams. Then, the detector states 1 1)1,2) 
are given by the state of the undetected photon repre- 
sented by diagram (2a) + (2c), on the one hand, and 
(26) + (2d), Fig. 3, on the other one. Similarly, we may 
also include the elastic component of the photodetection 
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FIG. 7. Intensity of the inelastic detection signal (suitably 
normalized) , as a function of the angle 6 between detector and 
incident light, for detuning 5 = V. In (a), both the final fre- 
quency lud and the positions ri2 of the elastically and inelas- 
tically scattering atom are fixed (u)d = w a t and r± = 2ri2/3), 
(b) shows the average over u)d, and (c) the average over both 
ujd and ri2. To make the loss of coherence due to the average 
over the power spectrum visible, the arrows in (b) indicate the 
contrast for perfect coherence, 71,11 = 1, but same asymmetry 
h/hi, cf. Eq. (69). The average over the positions, apart 
from washing out the side maxima, leads to a further loss of 
coherence, due to exchange of the inelastically and elastically 
scattering atom (see main text). 

signal, where the undetected photon is described by 
the diagrams in Fig. 2. Furthermore, also the relation 
between degree of coherence, 712, and enhancement fac- 
tor a, Eq. (79), remains valid for the total signal (as it is 
always the case if I\ = I 2 ): 



7 (e. + in) 



C 4 

— = OL — 1 = , 

L 4 + s 



(81) 



where we have inserted Eqs. (61,62) of the previous sec- 
tion. 

The average over the positions ri and r 2 of the atoms, 
finally, does not affect the intensity observed at 9 = 0; it 
only reduces the side maxima, and determines the shape 
of the backscattering 'cone'. As an example, Fig. 7(c) 
shows the result for an angular average over r 2 — ri, at 
fixed distance 7-12, where the cone shape is described by 
the function sin(x)/x. 

Finally, we want to stress that there is no loss of co- 
herence associated with the inelastic scattering 'on its 
own', but only in connection with the frequency filter- 
ing induced by the elastic scattering event. This can 
be demonstrated as follows: let us imagine that the re- 
sponse of the second atom is frequency-independent, i.e. 



B(oj) — const in Eq. (41). Then, the amplitudes of two 
reversed paths, see Eqs. (64,65), are identical, the unde- 
tected photon does not carry any which-path informa- 
tion, and we recover the enhancement factor two, even 
in the presence of inelastic scattering. Such a situation 
can be realized, e.g., by choosing atoms with different 
linewidths T 2 ^> T\, such that atom 2 cannot resolve the 
spectrum emitted by atom 1. In this significant 
reduction of the enhancement factor is observed only if 
we increase the distance 7*12 between the atoms, such that 
the propagation in the vacuum becomes relevant. 



D. Conclusion 

In summary, we have presented the first calculation 
of coherent backscattering in the presence of saturation. 
For two distant atoms, with single scattering excluded, 
the slope of the backscattering enhancement factor as a 
function of the incident intensity so at s — equals 
— 1/4, independently of the value of the detuning. The 
reduction of the enhancement factor can be traced back 
to the following two random processes: firstly, the fre- 
quency of the photons may be changed by the inelas- 
tic scattering event, which may, secondly, occur either 
at the first or at the second atom. Both processes (the 
latter one only for nonzero detuning, see Eq. (74)) lead 
to a random phase shift between the doubly scattered 
light emitted by the first atom, on the one hand, and by 
the second atom, on the other one, resulting in a loss of 
coherence. Alternatively, the coherence loss can be ex- 
plained by regarding the undetected photon as a which- 
path detector: its final state contains information about 
whether the detected photon has been emitted by the 
first or second atom, thereby partially destroying coher- 
ence between those paths. 

Starting from the solution of our model, we can think 
of extending it to more general scenarios in two differ- 
ent directions, either increasing the number of photons, 
to reach higher values of the saturation parameter, or 
the number of atoms, to treat a disordered medium of 
atoms. Since the complexity of the scattering approach 
increases dramatically with the number of scattered par- 
ticles, it may be more promising to use other methods, 
such as the optical Bloch equations [17], in the case of 
high saturation. The opposite is true for a large num- 
ber of scatterers, where we can resort to known concepts 
from the theory of multiple scattering. An important 
question, which must be solved in order to interpret the 
results of the experiment [16], is how the average prop- 
agation of the two-photon state in the atomic medium 
affects the coherent backscattering signal. 
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APPENDIX A: 



TWO-ATOM SCATTERING 
MATRIX 



In this appendix, we calculate the scattering of two 
photons by two atoms. For this purpose, we use the 
following expansion of the evolution operator U(t ,t) = 
exp[-i(flo + V)(t-t )]: 



oo ft pi 

U(t ,t) = J2 / dh / dt 

n = Jt Jt! 



dtn 



Uoito^JVUoitut^V . . . VU Q (t n , t), (Al) 

where Uo(to, t) = CTcp[—iHo(t — to)] denotes the free evo- 
lution. With each interaction V, see Eq. (3), an atom 
may emit a photon, or absorb one of the two photons. 
The corresponding 'paths' connecting the initial and final 
two-photon state \i) = |k 1 e 1 ,k 2 £2) and |/) = |k 3 e3,k 4 e4) 
can be represented diagrammatically, see Fig. 8. 



Here, (Ia,b) describes the scattering of two photons by 
a single atom [18], and (IIa,b) and (IIIa,b) the scattering 
by two atoms. Let us first concentrate on (Ia,b) and 
(IIa,b), where the inelastic scattering event occurs before 
the elastic one. Note that in (Ha), we have not specified 
the order in which the photons are emitted or absorbed. 
What we mean by this is a sum over all possible orderings, 
as indicated in Fig. 9. As we will see below, however, the 
sum need not be explicitly evaluated. 




X 



^x 



+ etc. 



FIG. 9. Independent scattering of two photon by two dif- 
ferent atoms. This diagram appears as a building block in 
(Ha) and (Illb), Fig. 8. 

Furthermore, we have selected one of the interaction 
operators V in the expansion (Al), at which we split the 
diagrams into a right and left half, denoted by J7('' r ) in 
the following. According to Eq. (Al), we may write: 



Uu&(to, t) 



-f 

J t 



dhU^t^VU^t), (A2) 



(r), 



and similarly for the other three diagrams (la), (lb), and 
(lib). Obviously, the left half is identical in the one- and 
two-atom case I and II, respectively. In the right half, on 
the other hand, the two photons are always independent 
from each other, being scattered by different atoms (if 
at all). This means that the evolution operator is the 
product of the two single-photon evolution operators: 




3 4 



Ha) t 

r xjx 




Ilia) 



2\ 



Illb) 



X . \ 



YX 



M ^3 ^3 ~ d *4 

FIG. 8. Diagrams describing scattering of two photons by 
a single atom (Ia,b) and two atoms, (IIa,b) and (IIIa,b). The 
curly lines represent photons, and thin or thick lines an atom 
in the ground or excited state. In order to simplify the com- 
parison between (I) and (II), we split the diagrams into a right 
and a left half (see main text). 



(A3) 



and likewise for (la), (lb), and (lib). Note that the evo- 
lution of the first photon (1 — > 3) depends only on (I) or 
(II), and not on (a) or (b), and vice versa for the second 
photon. Thereby, if we want to compare the one- and 
two-atom case, we have to consider only the two single- 

(r 1) 

photon diagrams U l { Y , which arc illustrated in Fig. 10. 



I) 




k 3 £ 3 



FIG. 10. (I) Emission of a photon ^363) at time t\. (II) 
Photon emission and subsequent scattering by the second 
atom. 

The first one simply describes the emission of photon 
|k3£3) by an atom located at iq, followed by free evolu- 
tion: 



VU^\t u t) 



ig{e 1 ■ el)e 



-ik 3 -ri —iu> 3 {t—ti) 



(A4) 
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In the second case, the photon is scattered by the other 
atom. Here, we have to take the sum over its intermediate 
state. For the calculation, it is convenient to express 
the time evolution in terms of the corresponding Green's 
function: 

VU^\ tl ,t) = j^^-e-^-^G^iz), (A5) 

where the contour C + runs just above the real axis, i.e., 
z = x + ie, e > 0, from x = +00 to —00, and the Green's 
function of the above diagram (II) reads: 



-«5 3 (ei -e*)(e-c5)e 



— ik-ri+i(k— k 3 )-r 2 



(Z - LOk)(z - LU )(Z - UJ 3 ) 



(A6) 



In the continuous limit (L — > 00), the sum is replaced 
by an integral (J2k = (L/2n) 3 J dk). The result of the 
integral (A6), in leading order of 1/(^3^2), reads: 



3iTg(e 1 ■ A 12 ■ e *) e - lk3 - r2 e izri2 
4w 3 r 12 (z - w )(z - u> 3 ) 



(A7) 



Here, A12 denotes the projection onto the plane orthogo- 
nal to r2 — ri. Finally, in the contour integral (A5), only 
the pole at z = u 3 contributes (if t — t\ S> 1/r): 



VU^itut) = ig{e 1 ■ Ai 2 • e^e-^e" 



•iw 3 (t— ti) 



4w 3 ri 2 (w 3 - uj ) 



(A8) 



Comparing Eqs. (A4,A8), we see that the contribution 
to the two-atom scattering matrix represented by (IIa,b) 
is given by the one- atom matrix Si, times a correction 
of the geometrical phase and the polarization, times the 
photon exchange factor B(u) 3 ), see Eq. (41). 

(k 3 e3,k 4 e4|S'ii|kiei,k 2 e 2 ) = (k 3 e 3 , k 4 e 4 |Si|kiei, k 2 e 2 ) 

ik 3 . (ri -r 2 ) (£I^12^|) (AQ) 



xe" 



What remains is the contribution, where the elastic scat- 
tering occurs before the inelastic one, represented by dia- 
grams (IIIa,b) in Fig. 8. The calculation can be repeated 
in almost the same way as above, or simply by noting 
that (IIIa,b) is related to (IIa,b) through time reversal, 
and the result is: 

(k 3 e 3 ,k 4 e 4 |S'iii|kiei,k2e2) = (k 3 e 3 , k 4 e 4 |S'i|kiei, k 2 e 2 ) 



xe 



-^•^-^) (ei ; Aia ;^ ) B(o; 1 ). (A10) 



Here, the photon exchange factor B(w) is evaluated at 
the frequency of the initial photon. The total scattering 
matrix is now readily obtained by adding Sn and Sm, 
and also including the diagrams where the two atoms 
and/or the two photons are exchanged. 



APPENDIX B: THE SCALAR CASE 

In this appendix, we calculate the photodetection sig- 
nal for scalar photons. Although they are not suited for 
coherent backscattering, since single scattering cannot 
be excluded, the solution will be useful for a future com- 
parison with the results obtained from the optical Bloch 
equations, which can be solved much more easily in the 
scalar case. 

As in the vectorial case, we consider contributions to 
the detection signal up to second order in l/(w_Lri 2 ). We 
neglect those terms whose order in 1/{u>lTyi) is changed 
by the angular average over ri 2 . 5 Furthermore, we con- 
sider only contributions which do not oscillate rapidly as 
a function of ri 2 , i.e. which survive an average over r\i 
over one wavelength. 

Firstly, since the two atoms may scatter independently 
from each other, we obtain two times the single-atom 
result, see Eqs. (34,35): 



=r ls (.s-2s 2 ), 

L (in,0) 



(Bl) 
(B2) 



which contributes to the background intensity L. Here 
we have to take into account that the lifetime T, and 
the prefactor 77, are different in the scalar and vectorial 
case, respectively. Instead of Eqs. (5,36), the following 
expressions hold for scalar photons: 



,3 



Vs 



2ne 
V 2dw L R 



(B3) 
(B4) 



Next, we consider the cases where one photon is ex- 
changed between the two atoms. These contribute to 
the detection signal in second order of l/(o>£ri 2 ). Con- 
cerning onc-photon scattering, only the diagrams (ld,e), 
Fig. 2, are relevant, and we obtain the same result as for 
the h || h channel in the vectorial case, see Eq. (46): 



£ (el,l) = C (el,l) = f?a | B |2 Si 

but with modified 'photon exchange factor' 



B 



compare Eq. (41). 



2LU L ri 2 (uJL - wo) 



(B5) 



(B6) 



5 These terms give the corrections of the average photon 
propagation induced by a disordered medium consisting of 
only a single atom. In the case of many disordered atoms, they 
are taken into account by renormalizing the single-photon 
propagation, in order to describe the mean free path and re- 
fractive index of the atomic medium. 
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The elastic contribution quadratic in s arises from in- 
terference of two-photon and one-photon scattering. Let 
us first look at diagram (2a). As before in the h \\ h chan- 
nel, it interferes with (la+lb) for the undetected photon, 
and (Id) or (le) for the detected photon, giving rise to 
— rj s \B\ 2 s 2 in background L and the cone C, respectively. 
Including single scattering, we obtain a new contribution: 
the detected photon may be singly scattered (lb), and 
the undetected photon either doubly scattered (le), or 
singly-scattered by the other atom (lc). Here, the state 
(le+lc), of the undetected photon exactly corresponds to 
the state (la+lb), in the previous case. Consequently, 
we obtain another term —ri s \B\ 2 s 2 in the background. 

With diagram (2b), the above considerations can be 
repeated in almost the same way. The difference from 
(2a) is only that the detected photon propagates in the 
opposite direction. Consequently, we obtain a new term 
— rj s \B\ 2 s 2 in the cone C, instead of the background L. 

Diagram (2e) is identical to diagram (2a), since we can- 
not distinguish between singly or doubly scattered pho- 
tons (open or full arrows in Figs. 2-4) in the scalar case. 
(2c), (2d), and (2f), finally, are obtained by exchanging 
the atoms. Adding all contributions mentioned above, 
we get: 



by subsequent scattering events (due to energy conserva- 
tion). Hence, the final result is given by Eqs. (B1-B10). 



L (el,2) = _1 0??S |,B|V 

C (cl ^ = -8 Vs \B\ 2 s 2 . 



(B7) 
(B8) 



As for the inelastic component, we only have to include 
the new diagrams (2e,f), which - as already mentioned 
above - are identical to (2a,b). Hence, the background 
contribution 2Ij, see Eq. (52), is multiplied by a factor 4, 
and the backscattering cone, Eq. (56), by a factor 2. We 
obtain: 

= 8h + 2J„ = P£ + f*) Vs\Bfs' (B9) 



C( in - 2 ) = 3 Vs \B\ 2 s 2 



(BIO) 



What we have not taken into account so far is interfer- 
ence between two diagrams where the undetected photon 
is emitted by different atoms. According to Eq. (49), the 
angular integral over the undetected photon then yields 
sin(o>ri2)/(wri 2 ). Hence, if one of the two diagrams con- 
tains a photon exchange, we obtain a contribution pro- 
portional to \B\ 2 . However, it can be shown that these 
contributions are exactly canceled by other contributions 
originating from the diagrams (2g,h), which also have 
been neglected so far. For example, the interference of 
(2g) with (lc) for the detected photon and (lc+lc) for 
the undetected one is canceled by the interference of (2j) 
with (lc) for the detected photon and (le) for the un- 
detected one. Similarly, the term |2g| 2 is canceled by 
the interference of (2g) with (2j). The underlying reason 
for all these cancellations is that what the undetected 
photon does after the inelastic scattering is irrelevant. 
We are only interested in its norm, which is not changed 
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